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(57) Abstract 

At least one parameter of at least one 
fluid in a pipe (12) is measured using a 
spatial array of acoustic pressure sensors 
(14, 16, 18) placed at predetermined ax- 
ial locations xl,x2,x3 along the pipe (12). 
The pressure sensors (14, 16, 18) provide 
acoustic pressure signals Pi(t), P2(t), P3(t) 
on lines (20, 22, 24) which are provided to 
signal processing logic (60) which deter- 
mines the speed of sound a m ix of the fluid 
(or mixture) in the pipe (12) using acoustic 
spatial array signal processing techniques 
with the direction of propagation of the 
acoustic signals along the longitudinal axis 
of the pipe (12). Numerous spatial array 
processing techniques may be employed to 
determined the speed of sound a m i x . The 
speed of sound a m j x is provided to logic 
(48) which calculates the percent compo- 
sition of the mixture, e.g., water fraction, 
or any other parameter of the mixture or 
fluid which is related to the sound speed 
a m u. The logic (60) may also determine 
the mach number (MX) of the fluid. The 
acoustic pressure signals Pi(t), P2(t), P3(t) 
measured are lower frequency (and longer 
wavelength) signals than those used for ul- 
trasonic flow meters, and thus is more tol- 
erant to inhomogeneities in the flow. No 
external source is required and thus may 
operate using passive listening. The invention will work with arbitrary sensor spacing and with as few as two sensors if certain information 
is known about the acoustic properties of the system. 




FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Airocnia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN ' 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


uz 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


zw 


Zimbabwe 


CI 


Cflte d'lvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


FT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







WO 00/00793 



PCT/US99/14589 



Fluid Parameter Measurement in Pipes Using Acoustic Pressures 

Cross References to Related Applications 

5 This application is an continuation-in-part of commonly owned co-pending 

US Patent Applications, Serial No., 09/105,534, entitled "Fluid Parameter 
Measurement in Pipes Using Acoustic Pressures", filed June 26, 1998, and contains 
subject matter related to that disclosed in commonly owned co-pending US Patent 
Applications: Serial No. (CiDRA Docket No. CC-0187), entitled "Measurement of 

10 Propagating Acoustic Waves in Compliant Pipes", filed contemporaneously herewith, 
Serial No. (CiDRA Docket No. CC-0194), entitled "Displacement Based Pressure 
Sensor Measuring Unsteady Pressure in a Pipe", and Serial No. (CiDRA Docket No. 
CC-0102A), entitled "Non-Intrusive Fiber Optic Pressure Sensor for Measuring 
Unsteady Pressures within a Pipe", filed contemporaneously herewith, all of which 

15 are incorporated herein by reference. 

Technical Field 

This invention relates to fluid parameter measurement in pipes and more 
particularly to measuring speed of sound and parameters related thereto of fluids in 
20 pipes using acoustic pressures. 

Background Art 

It is known that the speed of sound a^x of fluids in pipes may be used to 
determine various parameters of the fluid, such as is described in US Patent No. 

25 4,080,837, entitled "Sonic Measurement of Flow Rate and Water Content of Oil- 
Water Streams", to Alexander et al, US Patent No. 5,1 15,670, entitled "Measurement 
of Fluid Properties of Two-Phase Fluids Using an Ultrasonic Meter", to Shen, and US 
Patent 4,1 14,439, entitled "Apparatus for Ultrasonically Measuring Physical 
Parameters of Flowing Media", to Fick. Such techniques have a pair of acoustic 

30 transmitters/receivers (transceivers) which generate a sound signal and measure the 
time it takes for the sound signal to travel between the transceivers. This is also 
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known as a "sing-around" or "transit time" method. However, such techniques 
require precise control of the acoustic source and are costly and/or complex to 
implement in electronics. 

Also, these techniques use ultrasonic acoustic signals as the sound signal 
5 measured, which are high frequency, short wavelength signals (i.e., wavelengths that 
are short compared to the diameter of the pipe). Typical ultrasonic devices operate 
near 200k Hz, which corresponds to a wavelength of about 0.3 inches in water. In 
general, to allow for signal propagation through the fluid in an unimpeded and thus 
interpretable manner, the fluid should be homogeneous down to length scales of 

10 several times smaller than the acoustic signal wavelength. Thus, the criteria for 

homogeneity of the fluid becomes increasingly more strict with shorter wavelength 
signals. Consequently, inhomogeneities in the fluid, such as bubbles, gas, dirt, sand, 
slugs, stratification, globules of liquid, and the like, will reflect or scatter the 
transmitted ultrasonic signal. Such reflection and scattering inhibit the ability of the 

1 5 instrument to determine the propagation velocity. For this reason, the application of 
ultrasonic flowmeters have been limited primarily to well mixed flows. 

Summary of the Invention 

Objects of the present invention include provision of a system for measuring 
20 the speed of sound of fluids in pipes. 

According to the present invention, an apparatus for measuring at least one 
parameter of a mixture of at least one fluid in a pipe, comprising a spatial array of at 
least two pressure sensors, disposed at different axial locations along the pipe, and 
each measuring an acoustic pressure within the pipe at a corresponding axial location, 
25 each of said sensors providing an acoustic pressure signal indicative of the acoustic 
pressure within the pipe at said axial location of a corresponding one of said sensors; 
and a signal processor, responsive to said pressure signals, which provides a signal 
indicative of a speed of sound of the mixture in the pipe. 

According further to the present invention, the signal processor comprises 
30 logic which calculates a speed at which sound propagates along the spatial array. 
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According farther to the present invention, the signal processor comprises 
logic which calculates a frequency domain representation of (or frequency based 
signal for) each of the acoustic pressures signals. According still further to the 
present invention, the signal processor comprises logic which calculates a ratio of two 
5 of the frequency signals. In still further accord to the present invention, the sensors 
comprise at least three sensors. 

According still further to the present invention, the pressure sensors are fiber 
optic Bragg grating-based pressure sensors. Still further accord to the present 
invention, at least one of the pressure sensors measures an circumferential-averaged 

10 pressure at a given axial location of the sensor. Further according to the present 
invention, at least one of the pressure sensors measures pressure at more than one 
point around a circumference of the pipe at a given axial location of the sensor. 

The present invention provides a significant improvement over the prior art by 
providing a measurement of the speed of sound a mix of a mixture of one or more fluids 

15 within a pipe (where a fluid is defined as a liquid or a gas) by using an axial array of 
acoustic (or ac, dynamic, unsteady, or time varying) pressure measurements along the 
pipe. An explicit acoustic noise source is not required, as the background acoustic 
noises within the pipe (or fluid therein) will likely provide sufficient excitation to 
enable characterization of the speed of sound of the mixture by merely passive 

20 acoustic listening. 

The invention works with acoustic signals having lower frequencies (and thus 
longer wavelengths) than those used for ultrasonic meters, such as below about 20k 
Hz (depending on pipe diameter). As such, the invention is more tolerant to the 
introduction of gas, sand, slugs, or other inhomogeneities in the flow. 

25 The invention will work with arbitrary sensor spacing and arbitrary flow Mach 

numbers Mx; however, if the sensors are equally spaced and the axial velocity of the 
flow is small and therefore negligible compared to the speed of sound in the mixture 
(i.e., Mach number of the mixture Mx is small compared to one), the speed of sound 
a m j X may be determined as an explicit function of the frequency domain representation 

30 (frequency based signal) for the acoustic pressure signals at a given evaluation 
frequency go. 
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Since the speed of sound is an intrinsic property of mixtures, the present 
invention can be used to measure any parameter (or characteristic) of any mixture of 
one or more fluids in a pipe in which such parameter is related to the speed of sound 
of the mixture a,™, e.g., fluid fraction, temperature, salinity, sand particles, slugs, pipe 
5 properties, etc. or any other parameter of the mixture that is related to the speed of 
sound of the mixture. For example, the present invention may be used to measure 
fluid volume fractions (or composition or cut or content) of a mixture of any number 
of fluids in which the speed of sound of the mixture a^ is related to (or is 
substantially determined by), the volume fractions of two constituents of the mixture, 
10 e.g., oil/water, oil/gas, water/gas. Also, the present invention can be used to measure 
the speed of sound of any mixture and can then be used in combination with other 
known quantities to derive phase content of mixtures with multiple (more than two) 
constituents. 

The present invention allows the speed of sound to be determined in a pipe 

15 independent of pipe orientation, i.e., vertical, horizontal, or any orientation 

therebetween. Also, the invention does not require any disruption to the flow within 
the pipe (e.g., an orifice or venturi). Further, the invention uses ac (or unsteady or 
dynamic) pressure measurements as opposed to static (dc) pressure measurements and 
is therefore less sensitive to static shifts (or errors) in sensing. Furthermore, if harsh 

20 environment fiber optic pressure sensors are used to obtain the pressure 

measurements, such sensors eliminate the need for any electronic components down- 
hole, thereby improving reliability of the measurement. 

Also, a strain gauge (optical, electrical, etc.) that measures hoop strain on the 
pipe may be used to measure the ac pressure. Fiber optic wrapped sensors may be 

25 used as optical strain gauges to provide circumferentially-averaged pressure. Thus, 
the present invention provides non-intrusive measurements of the speed of sound (and 
other corresponding parameters), which enables real time monitoring and 
optimization for oil and gas exploration and production, or for other applications. 
The foregoing and other objects, features and advantages of the present 

30 invention will become more apparent in light of the following detailed description of 
exemplary embodiments thereof. 
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Brief Description of the Drawings 

Fig. 1 is a schematic block diagram of a fluid parameter measurement system, 
in accordance with the present invention. 
5 Fig. 2 is a graph of the speed of sound of a mixture versus the percent water 

volume fraction for an oil/water mixture, in accordance with the present invention. 

Fig. 3 is a transmission matrix model for the acoustics of an example pipe 
having 9 sections and a radiation impedance in accordance with the present 
invention. 

10 Fig. 4, illustrations (a)-(c), are graphs of axial values for p m i x , a m i x , h wat cr 

properties of a mixture for the segments of the pipe of Fig. 3, in accordance with the 
present invention. 

Fig. 5 is a graph of magnitude and phase versus frequency for a ratio of two 
pressures P1/P2, for radiation impedance of 1.0, water fraction of 50%, and axial 
1 5 properties of Fig. 4, in accordance with the present invention. 

Fig. 6 is a graph of magnitude and phase versus frequency for a ratio of two 
pressures P1/P3, for radiation impedance of 1.0, water fraction of 50%, and axial 
properties of Fig. 4, in accordance with the present invention. 

Fig. 7 is a graph of the magnitude of the speed of sound estimate versus an 
20 error term over a range of frequencies, using the frequency responses of Figs. 5,6, in 
accordance with the present invention. 

Fig. 8 is a graph of magnitude and phase versus frequency for a ratio of two 
pressures P1/P2, for radiation impedance of 0.5, water fraction of 50%, and constant 
axial properties of the mixture, in accordance with the present invention. 
25 Fig. 9 is a graph of magnitude and phase versus frequency for a ratio of two 

pressures P1/P3, for radiation impedance of 0.5, water fraction of 50%, and constant 
axial properties of the mixture, in accordance with the present invention. 

Fig. 10 is a graph of magnitude and phase versus frequency for a ratio of two 
pressures P1/P2, for radiation impedance of 0.5, water fraction of 5%, and constant 
30 axial properties of the mixture, in accordance with the present invention. 
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Fig. 1 1 is a graph of magnitude and phase versus frequency for a ratio of two 
pressures P1/P3, for radiation impedance of 0.5, water fraction of 5%, and constant 
axial properties of the mixture, in accordance with the present invention. 

Fig. 12 is a graph of the magnitude of the speed of sound estimate versus an 
5 error term over a range of frequencies, using the frequency response for two different 
percent water fractions, of Figs. 8-1 1, in accordance with the present invention 

Fig. 13 is a contour plot of speed of sound versus axial Mach versus an error 
term, for 5% water fraction, Mach number of 0.05, at 25 Hz, in accordance with the 
present invention. 

10 Fig. 14 is a contour plot of speed of sound versus axial Mach versus an error 

term, for 50% water fraction, Mach number of 0.05, at 25 Hz, in accordance with the 
present invention. 

Fig. 15 is a portion of a logic flow diagram for logic of Fig. 1, in accordance 
with the present invention. 
1 5 Fig. 1 6 is a continuation of the logic flow diagram of Fig. 1 5, in accordance 

with the present invention. 

Fig. 1 7 is a schematic block diagram of a fluid parameter measurement 
system, in an oil or gas well application, using fiber optic sensors, in accordance with 
the present invention. 

20 Fig. 1 8 is a plot of speed of sound against wall thickness of a pipe for a rigid 

and a non-rigid pipe, in accordance with the present invention. 

Fig. 19 is a cross-sectional view of a pipe, showing a plurality of sensors 
around the circumference of the pipe, in accordance with the present invention. 

Fig. 20 is a side view of a pipe having an isolating sleeve around the sensing 
25 region of the pipe, in accordance with the present invention. 

Fig. 21 is an end view of a pipe showing pressure inside and outside the pipe, 
in accordance with the present invention. 

Fig. 22 is a side view of a pipe having optical fiber wrapped around the pipe at 
each unsteady pressure measurement location and a pair of Bragg gratings around 
30 each optical wrap, in accordance with the present invention. 
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Fig. 23 is a side view of a pipe having optical fiber wrapped around the pipe at 
each unsteady pressure measurement location with a single Bragg grating between 
each pair of optical wraps, in accordance with the present invention. 

Fig. 24 is a side view of a pipe having optical fiber wrapped around the pipe at 
5 each unsteady pressure measurement location without Bragg gratings around each of 
the wraps, in accordance with the present invention. 

Fig. 25 is an alternative geometry of an optical wrap of Figs. 21,22, of a 
radiator tube geometry, in accordance with the present invention. 

Fig. 26 is an alternative geometry of an optical wrap of Figs. 21,22, of a race 
10 track geometry, in accordance with the present invention. 

Fig. 27 is a side view of a pipe having a pair of gratings at each axial sensing 
location, in accordance with the present invention. 

Fig. 28 is a side view of a pipe having a single grating at each axial sensing 
location, in accordance with the present invention. 
15 Fig. 29 is a top view of three alternative strain gauges, in accordance with the 

present invention. 

Fig. 30 is a side view of a pipe having three axially spaced strain gauges 
attached thereto, in accordance with the present invention. 

Fig. 31 is an end view of a pipe having three unsteady pressure sensors spaced 
20 apart from each other within the pipe, in accordance with the present invention. 

Fig. 32 is a side view of a pipe having three unsteady pressure sensors spaced 
axially within the pipe, in accordance with the present invention. 

Fig. 33 is a side view of a pipe having three unsteady pressure sensors axially 
and radially spaced within the pipe, in accordance with the present invention. 
25 Fig. 34 is a side view of a pipe having an inner tube with axially distributed 

optical fiber wraps for unsteady pressure sensors, in accordance with the present 
invention. 

Fig. 35 is a side view of a pipe having an inner tube with axially distributed 
unsteady pressure sensors located along the tube, in accordance with the present 
30 invention. 
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Fig. 36 is a side view of a pipe having an inner tube with three axially 
distributed hydrophones located within the tube, in accordance with the present 
invention. 

Fig. 37 is a diagram showing the propagation of acoustic waves from a single 
source in two dimensional space onto a spatial array, in accordance with the present 
invention. 

Fig. 38 is a side view of a pipe having left and right travelling acoustic waves 
propagating along the pipe, in accordance with the present invention. 

Fig. 39 is a diagram showing the propagation of acoustic waves from two 
sources in two dimensional space onto a spatial array, in accordance with the present 
invention. 

Fig. 40 is a schematic block diagram of an alternative embodiment of a fluid 
parameter measurement system, in accordance with the present invention. 

Fig. 41 is a graph of speed of sound versus water cut, in accordance with the 
15 present invention. 

Best Mode for Carrying Out the Invention 

Referring to Fig. 1, a pipe (or conduit) 12 has three acoustic pressure sensors 
14,16,18, located at three locations x,,x 2 ,x 3 along the pipe 12. The pressure may be 
measured through holes in the pipe 12 ported to external pressure sensors or by other 
techniques discussed hereinafter. The pressure sensors 14,16,18 provide pressure 
time-varying signals P,(t),P 2 (t),P 3 (t) on lines 20,22,24, to known Fast Fourier 
Transform (FFT) logics 26,28,30, respectively. The FFT logics 26,28,30 calculate the 
Fourier transform of the time-based input signals P,(t),P 2 (t), P 3 (t) and provide 
complex frequency domain (or frequency based) signals P,(co),P 2 (a>),P 3 (a)) on lines 
32,34,36 indicative of the frequency content of the input signals. Instead of FFT's, 
any other technique for obtaining the frequency domain characteristics of the signals 
Pi(t)^ 2 (t),P 3 (t), may be used. For example, the cross-spectral density and the power 
spectral density may be used to form a frequency domain transfer functions (or 
30 frequency response or ratios) discussed hereinafter. 
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Also, some or all of the functions within the logic 60 may be implemented in 
software (using a microprocessor or computer) and/or firmware, or may be 
implemented using analog and/or digital hardware, having sufficient memory, 
interfaces, and capacity to perform the functions described herein. 
5 The frequency signals Pi(co),P 2 ((o),P3(a)) are fed to a^-Mx Calculation Logic 

40 which provides a signal on a line 46 indicative of the speed of sound of the mixture 
amix (discussed more hereinafter). The a m i X signal is provided to map (or equation) 
logic 48, which converts a m i X to a percent composition of the fluid and provides a 
%Comp signal on a line 50 indicative thereof (as discussed hereinafter). Also, if the 

10 Mach number Mx is not negligible and is desired to be known, the calculation logic 
40 may also provide a signal Mx on a line 59 indicative of the Mach number. Mx (as 
discussed hereinafter). 

More specifically, for planar one-dimensional acoustic waves in a 
homogenous mixture, it is known that the acoustic pressure field P(x,t) at a location x 

15 along a pipe, where the wavelength A, of the acoustic waves to be measured is long 
compared to the diameter d of the pipe 12 (i.e., X/d »1), may be expressed as a 
superposition of a right traveling wave and a left traveling wave, as follows: 

P(x, t) = (Ae - ik ' x + Be + ik ' x }> im Eq. 1 

where A,B are the frequency-based complex amplitudes of the right and left traveling 
20 waves, respectively, x is the pressure measurement location along a pipe, co is 

frequency (in rad/sec, where co=27if), and k r ,ki are wave numbers for the right and left 
travelling waves, respectively, which are defined as: 



f \ 
CD 



nix J 



1 + M, 



and hi = 



' CO ^ 



\ a mix J 



Eq.2 



where a m i X is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), 
25 and M x is the axial Mach number of the flow of the mixture within the pipe, where: 

M x =^ Eq.3 

a mix 

where Vmix is the axial velocity of the mixture. For non-homogenous mixtures, the 
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axial Mach number represents the average velocity of the mixture and the low 
frequency acoustic field description remains substantially unaltered. 

The frequency domain representation P(x,oa) of the time-based acoustic 
pressure field P(x,t) within a pipe, is the coefficient of the e itat term of Eq. 1, as 
5 follows: 

P(x,a>) = Ae- ik r x +Be +ik ' x . Eq. 4 

Referring to Fig. 1, we have found that using Eq. 4 for P(x,a>) at three axially 
distributed pressure measurement locations x,,x 2> x 3 along the pipe 12 leads to an 
equation for a™ as a function of the ratio of frequency based pressure measurements, 
which allows the coefficients A,B to be eliminated. For optimal results, A and B are ' 
substantially constant over the measurement time and substantially no sound (or 
acoustic energy) is created or destroyed in the measurement section. The acoustic 
excitation enters the test section only through the ends of the test section 51 and, thus, 
the speed of sound within the test section 51 can be measured independent of the 
acoustic environment outside of the test section. In particular, the frequency domain 
pressure measurements P,(»)A(»)A(o») at the three locations x,,x 2 ,x 3 , respectively, 
along the pipe 12 using Eq. 1 for right and left traveling waves are as follows: 

Pi(co) =P(x = x t ,co) = Ae'** + Be**' x > Eq. 5 

P 2 (co) =P(x = x 2 ,co) = + Be +ik,x 2 Eq 6 

P3(co) =P(x = x 3 ,o)) = Ae** + Be + *"> Eq. 7 

where, for a given frequency, A and B are arbitrary constants describing the acoustic 
field between the sensors 14,16,18. Forming the ratio of P,((o)/P 2 (cd) from Eqns. 6,7, 
and solving for B/A, gives the following expression: 



R = - = 
A 



B_ [P 2 (o>) 



Eq.8 



P 2 (a>) 

25 where R is defined as the reflection coefficient. 

Forming the ratio of P,(co)/P3(g>) from Eqs. 5 and 7 and solving for zero gives: 
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+ Re 



ik,x, 



= 0 



Eq.9 



where R=B/A is defined by Eq. 8 and kr and kl are related to a m i X as defined by Eq. 2. 
Eq. 9 may be solved numerically, for example, by defining an "error" or residual term 
as the magnitude of the left side of Eq.9, and iterating to minimize the error term. 



mag 



+ Re 



P x {a>) 



= Error 



Eq. 10 



For many applications in the oil industry, the axial velocity of the flow in the 
pipe is small compared to the speed of sound in the mixture (i.e., the axial Mach 
number M x is small compared to one). For example, the axial velocity of the oil V 0 ii in 
a typical oil well is about 10 ft/sec and the speed of sound of oil aoii is about 4,000 
ft/sec. Thus, the Mach number Mx of a pure oil mixture is 0.0025 (V 0 ii/aoii=l 0/4,000), 
and Eq. 2 reduces to approximately: 

(O 



Eq. 11 



and the distinction between the wave numbers for the right and left traveling waves is 
eliminated. In that case (where Mx is negligible), since all of the variables in Eq. 10 
are known except for a m j X , the value for a m j X can be iteratively determined by 
evaluating the error term at a given frequency co and varying a,™ until the error term 
goes to zero. The value of amix at which the magnitude of the error term equals zero 
(or is a minimum), corresponds to the correct value of the speed of sound of the 
mixture a m j X . As Eq. 10 is a function of frequency co, the speed of sound a m i X at which 
the error goes to zero is the same for each frequency co evaluated (discussed more 
hereinafter). However, in practice, there may be some variation over certain 
frequencies due to other effects, e.g., pipe modes, non-acoustical pressure 
perturbation, discretization errors, etc., which may be filtered, windowed, averaged, 
etc. if desired (discussed more hereinafter). Furthermore, since each frequency is an 
independent measurement of the same parameter, the multiple measurements may be 
weighted averaged or filtered to provide a single more robust measurement of the 
speed of sound. 
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One example of how the speed of sound of the mixture a™ in the pipe 12 may 
be used is to determine the volume fraction of the mixture. In particular, the speed of 
sound of a mixture a mix of two fluids (where a fluid is defined herein as a liquid or a 
gas) in a pipe is in general related to the volume fraction of the two fluids. This 
relationship may be determined experimentally or analytically. For example, the 
speed of sound of a mixture may be expressed as follows: 



where a,,a 2 are the known speeds of sound, p,,p 2 are the known densities, and h,,h 2 
are the volume fractions of the two respective fluids, a mix is the speed of sound of the 
mixture, and the densities pi,p 2 of the two fluids are within about an order of 
magnitude (10:1) of each other. Other expressions relating the phase fraction to speed 
of sound may be used, being derived experimentally, analytically, or computationally. 

Referring to Fig. 2, where the fluid is an oil/water mixture, a curve 10 shows 
the speed of sound of the mixture a mix plotted as a function of water volume fraction 
using Eq. 12. For this illustrative example, the values used for density (p) and speed 
of sound (a) of oil and water are as follows: 

Density (p): p water = 1,000 kg/m 3 ; p oil = 700 kg/m 3 
Speed of sound (a): a water = 5,000 ft/sec; a*, = 4,000 ft/sec. 
The subscripts 1,2 of Eq. 12 assigned to the parameters for each fluid is arbitrary 
provided the notation used is consistent. Thus, if the speed of sound of the mixture 
amix is measured, the oil/water fraction may be determined. 

Referring to Fig. 3, to illustrate the concept by example, a transmission matrix 
model for the acoustics of an example pipe having 9 sections (or elements or 
segments) 1-9, an acoustic source 64, a radiation (or transmission) impedance Cm 
(£rad ■ = P/PmixamixH-mix) where p mix is an acoustic perturbation; Mx=0, and where the 
pressures Pi,P 2 ,P 3 are measured across test sections 5-6 and 6-7. For this example, 
each element is 1 meter long. 



a, 



mix 




Eq. 12 
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Depending on the application, an explicit acoustic noise source may or may 
not be required, as the background acoustic noises within the pipe may provide 
sufficient excitation to enable a speed of sound measurement from existing ambient 
acoustic pressures. In an oil or gas well application, if the background acoustic noises 
5 are not sufficient, an acoustic noise source (not shown) may be placed at the surface 
of the well or within the well, provided the source is acoustically coupled to the test 
section 51 over which the speed of sound is measured. 

Referring to Fig. 4, illustrations (a)-(c), an example of the axial properties of 
the mixture in the segments 1-9 of the pipe 12 is shown. The volume fraction of water 
10 h, the speed of sound of the mixture a mix , and the density of the mixture p m j x vary over 
the length of the pipe 12 and the test segments 5,6 (from 4-6 meters) between the 
pressure measurements P1-P3 have constant properties. In particular, the values for 
Pmix, a m j X , h wa ter for sections 1-9, respectively, are shown graphically in Fig. 4 and are 
as follows: 

15 h water = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9; 

Pmix = 730, 760, 790, 820, 850, 850, 910, 940, 970 (kg/m 3 ); 
amix =4053, 4111, 4177, 4251, 4334, 4334, 4539, 4667, 4818 (ft/sec); 
Referring to Figs. 5,6, the magnitude and phase of the ratio of the frequency 
based pressure signals Pi(co)/P 2 (co) and Pi(co)/P3(a>) is shown for the model of Fig. 3 
20 with the properties of Fig. 4 with 50% water in the test section and a radiation 
impedance of (^rad =1 .0 corresponding to an infinitely long pipe with constant 
properties of p m i X and a m i X for section 9 and beyond. 

Referring to Fig. 7, the error term of Eq. 10 using the frequency responses of 
Figs. 5,6, is a family of curves, one curve for each frequency co, where the value of 
25 the error is evaluated for values of a m j X varied from a wa t CT (5,000 ft/sec) to a^i (4,000 
ft/sec) at each frequency and the frequency is varied from 5 to 200 Hz in 5 Hz 
increments. Other frequencies may be used if desired. The speed of sound a m ix where 
the error goes to zero (or is minimized) is the same for each frequency co evaluated. 
In this case, the error is minimized at a point 70 when a m i X is 4335 ft/sec. From Fig. 2, 
30 for an oil/water mixture, an a mix of 4335 ft/sec corresponds to a 50% water volume 
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ratio in the test section which matches the water fraction of the model. Also, the 
sensitivity of a change in a™ to a change in error varies based on the evaluation 
frequency. Thus, the performance may be optimized by evaluating a™ at specific 
low sensitivity frequencies, such frequencies to be determined depending on the 
5 specific application and configuration. 

Referring to Figs. 8,9, for an radiation impedance £ad = 0.5, the magnitude 
and phase of the frequency responses (i.e., the ratio of frequency based pressure 
signals) P,(a>yp 2 (a>) and U<o)/P 3 (a>) is shown for the model of Fig. 3 with constant 
properties across all sections 1-9 of 50% water fraction (h=0.5), density of mixture 
10 Pmix= 850 kg/m 3 , and speed of sound of mixture a mix = 4334 ft/ sec . 

Referring to Fig. 12, for a 50% water fraction, the magnitude of the error term 
of Eq. 10 using the frequency responses of Figs. 8,9, is a family of curves, one curve 
for each frequency (o, where the value of a.* is varied from w (5,000 ft/sec) to a* 
(4,000 ft/sec) at each frequency and is shown at four frequencies 50,100,150,200 Hz. 
15 As discussed hereinbefore, the speed of sound a mix where the error goes to zero (or is 
minimized) is the same for each frequency co evaluated. In this case, the error is 
minimized at a point 72 where a™ = 4334 ft/sec, which matches the value of a™ 
shown in Fig. 7 for the same water fraction and different ^rad. From Fig. 2 (or Eq. 2), 
for an oil/water mixture, an a mix of 4334 ft/sec corresponds to a 50% water volume 
ratio in the test section which corresponds to the water fraction of the model. This 
shows that the invention will accurately determine a mix independent of the acoustic 
properties of the mixture outside the test sections and/or the termination impedances. 

Referring to Figs. 10,1 1, the magnitude and phase of the frequency responses 
(i.e., the ratio of the frequency based pressure signals) P,(co)/P 2 (a)) and P,(o))/P 3 (g)) is 
shown for the model of Fig. 3 with constant properties across all sections 1-9 of 5% 
water fraction (h=0.05), density of mixture Pmix = 715 kg/m 3 , and speed of sound of 
mixture a^ 4026 ft/sec, and a radiation impedance £ad = 0.5. 

Referring to Fig. 12, for a 5% water fraction, the magnitude of the error term 
of Eq. 10 using the frequency responses of Figs. 10,1 1, is a family of dashed curves, 
one curve for each frequency o>, where the value of a^ is varied from awater (5,000 
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ft/sec) to aoii (4,000 ft/sec) at each frequency and is shown at four frequencies 
50,100,150,200 Hz. As discussed hereinbefore, the speed of sound a m j X where the 
error goes to zero (or is minimized) is the same for each frequency co evaluated. In 
this case, the error is minimized at a point 74 when a m j x = 4026 ft/sec. From Fig. 1 (or 
Eq. 1), for an oil/water mixture, an a m j X of 4026 ft/sec corresponds to a 5% water 
volume ratio in the test section which corresponds to the water fraction of the model 
and, thus, verifies the results of the model. 

Referring to Fig. 12, for both 5% and 50% water fraction, the sensitivity of a 
change in a,™ to a change in error varies based on the evaluation frequency. In 
particular, for this example, of the four frequencies shown, the error approaches zero 
with the largest slope (AError/Aa mix ) for the 200 Hz curve, thereby making it easier to 
detect the value where the error goes to zero, and thus the value of a mix . Thus, 200 Hz 
would likely be a robust frequency to use to determine the speed of sound in this 
example. 

If the pressure sensors are equally spaced (i.e., xl-x2 = x3-x2 = Ax; or 
Axl=Ax2=Ax) and the axial Mach number Mx is small compared to one (and thus, 
kr=kl=k), Eq. 10 may be solved for k (and thus ami X ) in a closed- form solution as a 
function of the pressure frequency responses (or frequency based signal ratios) as 
follows: 



k = 



*mix L 



_1_ 

Ax 



zlog 



P n + P !3 P 12 +fc 2 2 +2P 13 P 12 2 +P ]3 2 P l2 2 -4P n 2 )' 2 



2P„ 



Solving for a mix , gives: 



co 



' 1 " 




/'log 


_Ax_ 



-4P 13 r 



2P, 



13 



Eq. 13 



Eq. 14 



where P 12 =Pi(co)/P2(co), Pn =Pi(a))/P 3 (a)), i is the square root of-1, and the result of 
the Log[] function is an imaginary number, yielding a real number for the speed of 
sound amix. 

The analytical solution to the Eq. 10 shown in Eqs. 13 and 14 is valid 
primarily for the frequencies for which the length of the test section 5 1 along the pipe 
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12 (i.e., x3-xl or 2Ax for equally spaced sensors) is shorter than the wavelength X of 
the acoustic waves to be measured. This restriction is due to multiple possible 
solutions to the Eq. 10. Alternative solutions to Eq. 10 for other frequency ranges may 
be derived using a variety of known techniques. 

An alternative closed form solution for a mix (in a trigonometric form) from the 
three pressure Eqs. 5-7, where the pressure sensors are equally spaced and Mx is 
negligible (i.e, kl=kr), is as follows. Forming the ratio [P,(o>) + P 3 (<d)]/P 2 (©) from 
Eqs. 5-7, gives the following expression: 

/>(6>) + P 3 (a>) _ Ae-*** +Be +ikx > + Ae'**' + Be +ikx > 

P 2 (a>) ~ Ae^+Be^ ^ 15 

For equally spaced sensors, xl=0,x2=Ax, x3=2Ax (xl=0 for convenience 
only), which gives: 

P^ + P^o) A + B + Ae~ 2iktx 



= +Be 

P 2 (<*>) " Ae-* 6 * +Be +ikAl 
Dividing the numerator and denominator by A, gives: 

Pi(co) + P 3 {a>) 1 + R + e- 2UUa + Re +2ik6x 



P 2 (*>) 



Eq. 16 



Eq. 17 



where R=B/A is defined by Eq. 8 with xl=0,x2=Ax, which gives: 



A 



1 - 


>>(<»)] 

_P 2 (a>)_ 




P x {a>) 

AM. 





Eq. 18 



Plugging R into Eq. 17, gives: 
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P,(fi>) + /» 3 (fl?) 

P 2 (a>) 



1- 



P 2 {a>) 



-ikAx 



_ x 2 

/>(*>) 



-1 



+ (l + e +2 ^) 



-ikAx . 

e + 
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P 2 (fi?) 
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P 2 (o>) 
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Eq. 19 



+ikAx 



Simplifying Eq. 19, gives: 



P l (a>)+P i {a>) _ 
P 2 (a>) 



L e «* _ij(e-^) + (i_|.^ |( € — ) 



Eq. 20 



Distributing terms and simplifying, gives: 
P^ + P^a) _-e 



-2i*Ax + g+2'*Ai 



P 2 (G>) 



— e + e 



Eq.21 



Using the relation between exponents and the sine function, gives: 

P x {g)) + P 3 (cd) _ 2/sin(2ftAx) _ 2sin(foc)cos(fcc) 
P 2 {(D) 2i sin(kx) sin(fcc) 

Simplifying and substituting k=co/a m i X5 gives: 

= 2cos(£Ax) = 2cos 



Eq. 22 



Eq. 23 



10 Eq. 23 is particularly useful due to its simple geometric form, from which a,™ 

can be easily interpreted. In particular, a m j X can be determined directly by inspection 
from a digital signal analyzer (or other similar instrument) set up to provide a display 
indicative of the left side of Eq. 23, which will be a cosine curve from which a mix may 
be readily obtained. For example, at the zero crossing of the cosine wave, Eq. 23 will 

15 be equal to zero and a mix will be equal to 2co AX/71. Alternatively, Eq. 23 may be used 
to determine amix using an iterative approach where a measured function is calculated 
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from the left side of Eq. 23 (using the measured pressures) and compared to a cosine 
curve of the right side of Eq. 23 where amix is varied until it substantially matches the 
measured function. Various other curve fitting, parameter identification, and/or 
minimum error or solution techniques may be used to determine the value of amix 
5 that provides the best fit to satisfy Eq. 23. 

Solving Eq. 23 for a mix , gives the following closed-form solution: 
a . - ate 



10 



15 



mix 

-I 



cos" 



^m±w] cos -.irw + w) Eq - 24 



V 



2P 2 (o>) 



2 



Referring to Fig. 41, a graph of speed of sound (a mix ) versus water cut is 
shown where a mix is calculated using Eq. 23 as described hereinbefore. Fig. 41 is for 
a Schedule 160 steel pipe having a 2 inch ID, AX=2ft even spacing between three 
axial sensing locations, each sensor being a piezo-electric ac pressure sensor, there 
being four evenly deferentially spaced sensors at each axial sensing location. The 
line 452 shows the theoretical value for water cut based on Eq.12 and Fig. 2 discussed 
hereinbefore, and the circles are the calculated values for a mix . 

Alternatively, Eq. 9 may be written in trigonometric forni for arbitrary spacing 
between the pressure sensors and where Mx is negligible (kl=kr), as follows: 



sinl ^-(* 3 -*,) Lp 32 Sin f JL (X2 _ Xi) 

\ mix J I a. 



sin 



\ a mix J 



~—(x 3 -x 2 )\ = 0 Eq25 

) 

where P 32 = P 3 (o))/P 2 ((d) and P 12 = P,'(co)/P 2 (co). 

Referring to Figs. 13,14, if Mach number Mx is not negligible and/or is 

20 desired to be calculated, the value of Mx and a mix where the error term of Eq. 10 is 
zero can be uniquely determined from Eq. 10 for a given water fraction. In particular, 
for a given %water fraction, there is a unique value indicated by points 90,92 for 5% 
and 50% water cut, respectively. Known software search algorithms may be used to 
vary ^ and Mx over predetermined ranges to find the value of Mx and a™ where 

25 the error = 0 (discussed more hereinafter). 

Referring to Fig. 15, the calculation logic 40 begins at a step 100 where P 12 is 
calculated as the ratio of P l( co)/P 2 (o», and a step 1 02 where P 13 is calculated as the 
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ratio of Pi(co)/P3(a)). Next a step 103 determines whether the Mach number Mx of the 
mixture is negligible (or whether it is desirable to calculate Mx). If Mx is negligible, 
a step 104 determines if the sensors 14,16,18 are equally spaced (i.e., xl-x2=x2- 
x3=Ax). If equally spaced sensors, steps 106 set initial values for g>= ool (e.g., 100 
5 Hz) and a counter n=l . Next, a step 108 calculates a m j x (n) from the closed form 

solution of Eq. 14. Then, a step 1 10 checks whether the logic 40 has calculated a,™ at 
a predetermined number of frequencies, e.g., 10. If n is not greater than 10, steps 
1 1 2,1 14, increments the counter n by one and increases the frequency co by a 
predetermined amount (e.g., 10 Hz) and the step 108 is repeated. If the logic 40 has 

10 calculated a mix at 10 frequencies, the result of the step 1 16 would be yes and the logic 
40 goes to a step 116 which determines an average value for a m j x using the values of 
a m ix(n) over the 10 frequencies, and the logic 40 exits. 

If the sensors are not equally spaced, steps 120 set xl,x2,x3 to the current 
pressure sensor spacing, and set initial values for co= col (e.g., 100 Hz) and the 

15 counter n=l. Next, a step 122 sets a m i x = a m i X -min (e.g., aojf=4000 ft/sec) and a step 124 
calculates the error term from Eq. 10. Then, a step 126 checks whether error =. 0. If 
the error does not equal zero, a mix is incremented by a predetermined amount and the 
logic 40 goes to a step 124. 

If the error=0 (or a minimum) in step 126, a step 130 sets a m i X (n) =ami X . Next, a 

20 step 132 checks whether n is greater than or equal to 10. If not, a step 134 increments 
n by one and increases the frequency co by a predetermined amount (e.g., 10 Hz). If n 
is greater than or equal to 10, a step 138 calculates an average value for a m j X over the 
10 frequencies. 

Referring to Fig. 1 6, if the Mach number Mx is not negligible, steps 200-204 
25 sets initial conditions: ©=g>1 (e.g., 100 Hz); Mx=Mx-min (e.g., 0); an^a^x-,™ (e.g., 
aoii^OOO ft/sec). Then, a step 206 calculates the error term of Eq. 10 at a step 202. 
Next, a step 208 checks whether the error=0 (or a minimum). If not, a step 210 checks 
whether a^a^x-max (e.g., awater^OOO ft/sec). 

If the result of step 210 is no, a step 212 increases a m i x by a predetermined 
30 amount (e.g., 1 ft/sec) and the logic goes to step 206. If the result of step 210 is yes, a 
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step 214 increases Mx by a predetermined amount (e.g., 1) and the logic goes to step 
204. 

When step 208 indicates error=0 (or a minimum), a step 216 sets ^(nK, 
and Mx(n)=Mx, and a step 218 checks whether the values of amix and Mx have been 
5 calculated at 10 different frequencies. If not, a step 220 increments the counter n by 
one and a step 222 increases the value of the frequency co by a predetermined amount 
(e.g., 10 Hz). If the values of a mix and Mx have been calculated at 10 different 
frequencies (i.e., n is equal to 10), a step 224 calculates a average values for a mix (n) 
and Mx(n) at the 10 different frequencies to calculate a mix and Mx. The value for 
amix above is similar to that shown in Figs. 13,14, discussed hereinbefore, where the 
final value of a™ are the points 90,92 where the error equals zero. 

Instead of calculating an average value for a™ in steps 1 16,138,24, a mix may 
be calculated by filtering or windowing a^n), from predetermined frequencies. The 
number of frequencies and the frequencies evaluated may be any desired number and 
values. Also, instead of calculating a™ and/or Mx at more than one frequency, it may 
be calculated at only one frequency. Further, the logic shown in Figs. 15,16 isone of 
many possible algorithms to calculate a™ using the teachings herein. 

Referring to Figs. 1 and 18, the compliance (or flexibility) of the pipe 12 (or 
conduit) in the sensing region may influence the accuracy or interpretation of the 
measured speed of sound a™ of the mixture in two primary ways. 

Regarding the first way, referring to Fig. 18, flexing of the pipe 12 in the 
sensing region reduces the measured speed of sound a™ from the sound in an 
unbounded domain. The sound speed in an unbounded domain (infinite media) is a 
property that is closely linked with the fluid properties. In particular, the influence of 
pipe wall thickness (or compliance of the pipe) on measured speed of sound due 
reduction in the speed of sound for a pipe having a 2 inch nominal diameter and 
having 100% water (p w =l,000 kg/m 3 ; a w =5,000 ft/sec) inside the pipe and a vacuum 
(or air) outside the pipe diameter, is shown. The speed of sound of water in an. 
infinitely rigid pipe (i.e., infinite modulus) is indicated by a flat curve 350, and the 
speed of sound of water in a steel pipe is indicated by a curve 352. A point 354 on the 
curve 352 indicates the value of the speed of sound of about 4768 ft/sec for a 
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Schedule 80 steel pipe. Accordingly, the thicker the pipe wall, the closer the speed of 
sound approaches the value of 5,000 ft/sec for an infinitely rigid pipe. 

The errors (or boundary effects) shown in Fig. 18 introduced into the 
measurements by a non-rigid (or compliant) pipe 12 can be calibrated and corrected 
5 for to accurately determine the speed of sound in the fluid in an unbounded media. 
Thus, in this case, while the system (pipe) does modify the propagation velocity, such 
velocity can be mapped to the propagation velocity in an infinite media in a 
predictable fashion. 

In particular, for fluids contained in a compliant pipe, the propagation velocity 
10 of compression waves is influenced by the structural properties of the pipe. For a fluid 
contained in the pipe 12 surrounded with a fluid of negligible acoustic impedance 
(pa), the propagation velocity is related to the infinite fluid domain speed of sound 
and the structural properties via the following relation: 

- = — + a whereas — Eq. 26 

p mix a measured p mix Cl mix Et 

1 5 where R= the pipe radius, t is the pipe wall thickness, p m i X is the density of the 

mixture (or fluid), amix is the actual speed of sound of the mixture, a mea sured is the 
measured speed of sound of the mixture contained in the pipe 12, and E is the 
Young's modulus for the pipe material. Eq. 26 holds primarily for frequencies where 
the wavelength of the acoustics is long (e.g., greater than about 2 to 1) compared to 

20 the diameter of the pipe and for frequencies which are low compared to the natural 
frequency of the breathing mode of the pipe. Eq. 26 also applies primarily to 
wavelengths which are long enough such that hoop stiffness dominates the radial 
deflections of the pipe. 

For Fig. 18, the curve 352 (for 100% water) would be one of a family of 

25 curves for various different oil/water mixtures. For Eq. 26, the terms may be defined 
in terms of the density of each constituent, and the volumetric phase fraction as 
follows: 

1 N <h N N 

= Z— 4" where: A»«=I>iA andj>,=l 



p mix a 2 mix M pflr M isl 
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where p, is the density of the i* constituent of a multi-component mixture, a, is 
the sound speed of the i th constituent of the mixture, +, is the volumetric phase fraction 
of the i* constituent of the mixture, and N is the number of components of the 
mixture. Knowing the pipe properties, the densities and the sound speed (in an infinite 
5 domain) of the individual constituents, and the measured sound speed of the mixture, 
Eq. 26 can be solved for a™,,. Thus, amix can be detennined for a compliant pipe. The 
calibration of the pipe can be derived from other equations or from a variety of other 
means, such as analytical, experimental, or computational. 

For certain types of pressure sensors, e.g., pipe strain sensors, accelerometers, 
1 0 velocity sensors or displacement sensors, discussed hereinafter, it may be desirable 
for the pipe 12 to exhibit a certain amount of pipe compliance. 

Alternatively, to minimize these error effects (and the need for the 
corresponding calibration) caused by pipe compliance, the axial test section 5 1 of the 
pipe 12 along where the sensors 14,16,18 are located may be made as rigid as 
15 possible. To achieve the desired rigidity, the thickness of the wall 53 of the test 

section 51 may be made to have a predetermined thickness, or the test section 51 may 
be made of a very rigid material, e.g., steel, titanium, Kevlar®, ceramic, or other 
material with a high modulus. 

Regarding the second way, if the pipe 12 is compliant and acoustically 
20 coupled to fluids and materials outside the pipe 12 in the sensing region, such as the 
annulus fluid, casing, rock formations, etc., the acoustic properties of these fluids and 
materials outside the pipe 12 diameter may influence the measured speed of sound. 
Because the acoustic properties of such fluids and materials are variable and 
unknown, their affect on measured speed of sound cannot be robustly corrected by 
25 calibration (nor mapped to the propagation velocity in an infinite media in a 
predictable fashion). 

Referring to Fig. 20, to alleviate this effect, an outer isolation sleeve 410 (or 
sheath, shell, housing, or cover) which is attached to the outer surface of pipe 12 over 
where the pressure sensors 14,16,18 are located on the pipe 12. The sleeve 410 forms 
30 a closed chamber 412 between the pipe 12 and the sleeve 41 0. We have found that 
when the chamber 412 is filled with a gas such as air, the acoustic energy in the pipe 
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is not acoustically coupled to fluids and materials outside the pipe 12 in the sensing 
region. As such, for a compliant pipe the speed of sound can be calibrated to the 
actual speed of sound in the fluid in the pipe 12 as discussed hereinbefore. The sleeve 
410 is similar to that Serial No. (Cidra Docket No. CC-0187), entitled "Measurement 
5 of Propagating Acoustic Waves in Compliant Pipes", filed contemporaneously 
herewith, which is incorporated herein by reference. 

Referring to Fig. 19, instead of single point pressure sensors 14,16,18, at the 
axial locations xl,x2,x3 along the pipe 12, two or more pressure sensors, e.g., four 
sensors 400-406, may be used around the circumference of the pipe 12 at each of the 

10 axial locations xl,x2,x3. The signals from the pressure sensors 400-406 around the 
circumference at a given axial location may be averaged to provide a cross-sectional 
(or circumference) averaged unsteady acoustic pressure measurement. Other numbers 
of acoustic pressure sensors and annular spacing may be used. Averaging multiple 
annular pressure sensors reduces noises from disturbances and pipe vibrations and 

15 other sources of noise not related to the one-dimensional acoustic pressure waves in 
the pipe 12, thereby creating a spatial array of pressure sensors to help characterize 
the one-dimensional sound field within the pipe 12. 

The pressure sensors 14,16,18 described herein may be any type of pressure 
sensor, capable of measuring the unsteady (or ac or dynamic ) pressures within a pipe, 

20 such as piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), 

accelerometers (or geophones), velocity measuring devices, displacement measuring 
devices, etc. If optical pressure sensors are used, the sensors 14-18 may be Bragg 
grating based pressure sensors, such as that described in copending US Patent 
Application, Serial No. 08/925,598, entitled " High Sensitivity Fiber Optic Pressure 

25 Sensor For Use In Harsh Environments", filed Sept. 8, 1997. Alternatively, the 

sensors 14-18 may be electrical or optical strain gages attached to or embedded in the 
outer or inner wall of the pipe which measure pipe wall strain, including microphones, 
hydrophones, or any other sensor capable of measuring the unsteady pressures within 
the pipe 12. In an embodiment of the present invention that utilizes fiber optics as the 

30 pressure sensors 14-18, they may be connected individually or may be multiplexed 
along one or more optical fibers using wavelength division multiplexing (WDM), 
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time division multiplexing (TDM), or any other optical multiplexing techniques 
(discussed more hereinafter). 

Referring to Fig. 21, if a strain gage is used as one or more of the pressure 
sensors 14-18, it may measure the unsteady (or dynamic or ac) pressure variations Pin 
5 inside the pipe 12 by measuring the elastic expansion and contraction, as represented 
by arrows 350, of the diameter (and thus the circumference as represented by arrows 
351) of the pipe 12. In general, the strain gages would measure the pipe wall 
deflection in any direction in response to unsteady pressure signals inside the pipe 12. 
The elastic expansion and contraction of pipe 12 is measured at the location of the 
10 strain gage as the internal pressure P in changes, and thus measures the local strain 
(axial strain, hoop strain or off axis strain), caused by deflections in the directions 
indicated by arrows 351, on the pipe 12. The amount of change in the circumference 
is variously determined by the hoop strength of the pipe 12, the internal pressure P^, 
the external pressure P out outside the pipe 12, the thickness T w of the pipe wall 352, 
and the rigidity or modulus of the pipe material. Thus, the thickness of the pipe wall 
352 and the pipe material in the sensor sections 51 (Fig. 1) may be set based on the 
desired sensitivity of sensors 14-18 and other factors and may be different from the 
wall thickness or material of the pipe 12 outside the sensing region 51. 

Still with reference to Fig. 21 and Fig. 1, if an accelerometer is used as one or 
more of the pressure sensors 14-18, it may measure the unsteady (or dynamic or ac) 
pressure variations P in inside the pipe 12 by measuring the acceleration of the surface 
of pipe 12 in a radial direction, as represented by arrows 350. The acceleration of the 
surface of pipe 12 is measured at the location of the accelerometer as the internal 
pressure P in changes and thus measures the local elastic dynamic radial response of 
the wall 352 of the pipe. The magnitude of the acceleration is variously determined by 
the hoop strength of the pipe 12, the internal pressure P in , the external pressure P out 
outside the pipe 12, the thickness T w of the pipe wall 352, and the rigidity or modulus 
of the pipe material. Thus, the thickness of the pipe wall 352 and the pipe material in 
the sensing section 51 (Fig. 1) may be set based on the desired sensitivity of sensors 
30 14-18 and other factors and may be different from the wall thickness or material of 
the pipe 12 outside the sensing region 14. Alternatively, the pressure sensors 14-18 
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may comprise a radial velocity or displacement measurement device capable of 
measuring the radial displacement characteristics of wall 352 of pipe 12 in response to 
pressure changes caused by unsteady pressure signals in the pipe 12. The 
accelerometer, velocity or displacement sensors may be similar to those described in 
5 commonly-owned copending US Patent Application, Serial No. (CiDRA Docket No. 
CC-0194), entitled "Displacement Based Pressure Sensor Measuring Unsteady 
Pressure in a Pipe", filed contemporaneously herewith and incorporated herein by 
reference. 

Referring to Figs. 22,23,24, if an optical strain gage is used, the ac pressure 

10 sensors 14-18 may be configured using an optical fiber 300 that is coiled or wrapped 
around and attached to the pipe 12 at each of the pressure sensor locations as 
indicated by the coils or wraps 302,304,306 for the pressures Pi,P2,P3> respectively. 
The fiber wraps 302-306 are wrapped around the pipe 12 such that the length of each 
of the fiber wraps 302-306 changes with changes in the pipe hoop strain in response 

15 to unsteady pressure variations within the pipe 12 and thus internal pipe pressure is 
measured at the respective axial location. Such fiber length changes are measured 
using known optical measurement techniques as discussed hereinafter. Each of the 
wraps measure substantially the circumferentially averaged pressure within the pipe 
12 at a corresponding axial location on the pipe 12. Also, the wraps provide axially 

20 averaged pressure over the axial length of a given wrap. While the structure of the 
pipe 12 provides some spatial filtering of short wavelength disturbances, we have 
found that the basic principle of operation of the invention remains substantially the 
same as that for the point sensors described hereinbefore. 

Referring to Fig. 22, for embodiments of the present invention where the 

25 wraps 302,304,306 are connected in series, pairs of Bragg gratings 

(310,312),(314,316), (318,320) may be located along the fiber 300 at opposite ends of 
each of the wraps 302,304,306, respectively. The grating pairs are used to multiplex 
the pressure signals Pi,P2,P3 to identify the individual wraps from optical return 
signals. The first pair of gratings 310,312 around the wrap 302 may have a common 

30 reflection wavelength X\ 9 and the second pair of gratings 314,316 around the wrap 
304 may have a common reflection wavelength A,2, but different from that of the first 
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pair of gratings 310,312. Similarly, the third pair of gratings 318,320 around the wrap 
306 have a common reflection wavelength k 3 , which is different from Jl,,Jl 2 . 

Referring to Fig. 23, instead of having a different pair of reflection 
wavelengths associated with each wrap, a series of Bragg gratings 360-366 with only 
5 one grating between each of the wraps 302-306 may be used each having a common 
reflection wavlength X] . 

Referring to Figs. 22 and 23 the wraps 302-306 with the gratings 310-320 
(Fig.22) or with the gratings 360-366 (Fig.23) may be configured in numerous known 
ways to precisely measure the fiber length or change in fiber length, such as an 

10 interferometric, Fabry Perot, time-of-flight, or other known arrangements. An 

example of a Fabry Perot technique is described in US Patent. No. 4,950,883 "Fiber 
Optic Sensor Arrangement Having Reflective Gratings Responsive to Particular 
Wavelengths", to Glenn. One example of time-of-flight (or Time-Division- 
Multiplexing; TDM) would be where an optical pulse having a wavelength is 

15 launched down the fiber 300 and a series of optical pulses are reflected back along the 
fiber 300. The length of each wrap can then be determined by the time delay between 
each return pulse. 

Alternatively, a portion or all of the fiber between the gratings (or including 
the gratings, or the entire fiber, if desired) may be doped with a rare earth dopant 

20 (such as erbium) to create a tunable fiber laser, such as is described in US Patent No. 
5,317,576, "Continuously Tunable Single Mode Rare-Earth Doped Laser 
Arrangement", to Ball et al or US Patent No. 5,513,913, "Active Multipoint Fiber 
Laser Sensor", to Ball et al, or US Patent No. 5,564,832, "Birefringent Active Fiber 
Laser Sensor", to Ball et al, which are incorporated herein by reference. 

25 While the gratings 310-320 are shown oriented axially with respect to pipe 12, 

in Figs. 22,23, they may be oriented along the pipe 12 axially, circumferentially, or in' 
any other orientations. Depending on the orientation, the grating may measure 
deformations in the pipe wall 352 with varying levels of sensitivity. If the grating 
reflection wavelength varies with internal pressure changes, such variation may be 

30 desired for certain configurations (e.g., fiber lasers) or may be compensated for in the 
optical instrumentation for other configurations, e.g., by allowing for a predetermined 
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range in reflection wavelength shift for each pair of gratings. Alternatively, instead of 
each of the wraps being connected in series, they may be connected in parallel, e.g., 
by using optical couplers (not shown) prior to each of the wraps, each coupled to the 
common fiber 300. 

5 Referring to Fig. 24, alternatively, the sensors 14-18 may also be formed as a 

purely interferometric sensor by wrapping the pipe 12 with the wraps 302-306 without 
using Bragg gratings where separate fibers 330,332,334 may be fed to the separate 
wraps 302,304,306, respectively. In this particular embodiment, known 
interferometric techniques may be used to determine the length or change in length of 

10 the fiber 10 around the pipe 12 due to pressure changes, such as Mach Zehnder or 

Michaelson Interferometric techniques, such as that described in US Patent 5,218,197, 
entitled "Method and Apparatus for the Non-invasive Measurement of Pressure Inside 
Pipes Using a Fiber Optic Interferometer Sensor" to Carroll. The inteferometric wraps 
may be multiplexed such as is described in Dandridge, et al, 'Tiber Optic Sensors for 

15 Navy Applications", IEEE, Feb. 1991, or Dandridge, et al, "Multiplexed 

Intereferometric Fiber Sensor Arrays", SPIE, Vol. 1586, 1991, ppl76-183. Other 
techniques to determine the change in fiber length may be used. Also, reference 
optical coils (not shown) may be used for certain interferometric approaches and may 
also be located on or around the pipe 12 but may be designed to be insensitive to 

20 pressure variations. 

Referring to Figs. 25 and 26, instead of the wraps 302-306 being optical fiber 
coils wrapped completely around the pipe 12, the wraps 302-306 may have alternative 
geometries, such as a "radiator coil" geometry (Fig 25) or a "race-track" geometry 
(Fig. 26), which are shown in a side view as if the pipe 12 is cut axially and laid flat. 

25 In this particular embodiment, the wraps 302-206 are not necessarily wrapped 360 
degrees around the pipe, but may be disposed over a predetermined portion of the 
circumference of the pipe 12, and have a length long enough to optically detect the 
changes to the pipe circumference. Other geometries for the wraps may be used if 
desired. Also, for any geometry of the wraps described herein, more than one layer of 

30 fiber may be used depending on the overall fiber length desired. The desired axial 
length of any particular wrap is set depending on the characteristics of the ac pressure 
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desired to be measured, for example the axial length of the pressure disturbance 
caused by a vortex to be measured. 

Referring to Figs. 27 and 28, embodiments of the present invention include 
configurations wherein instead of using the wraps 302-306, the fiber 300 may have 
5 shorter sections that are disposed around at least a portion of the circumference of the 
pipe 12 that can optically detect changes to the pipe circumference. It is further 
within the scope of the present invention that sensors may comprise an optical fiber 
300 disposed in a helical pattern (not shown) about pipe 12. As discussed herein 
above, the orientation of the strain sensing element will vary the sensitivity to 
10 deflections in pipe wall 352 deformations caused by unsteady pressure signals in the 
pipe 12. 

Referring to Fig. 27, in particular, the pairs of Bragg gratings (310,312), 
(314,316), (318,320) are located along the fiber 300 with sections 380-384 of the fiber 
300 between each of the grating pairs, respectively. In that case, known Fabry Perot 
15 interferometric, time-of-flight or fiber laser sensing techniques may be used to 
measure the strain in the pipe, in a manner similar to that described in the 
aforementioned references. 

Referring to Fig. 28, alternatively, individual gratings 370-374 maybe 
disposed on the pipe and used to sense the unsteady variations in strain in the pipe 12 
20 (and thus the unsteady pressure within the pipe) at the sensing locations. When a 
single grating is used per sensor, the grating reflection wavelength shift will be 
indicative of changes in pipe diameter and thus pressure. 

Any other technique or configuration for an optical strain gage may be used: 
The type of optical strain gage technique and optical signal analysis approach is not 
25 critical to the present invention, and the scope of the invention is not intended to be 
limited to any particular technique or approach. 

For any of the embodiments described herein, the pressure sensors, including 
electrical strain gages, optical fibers and/or gratings among others as described herein 
may be attached to the pipe by adhesive, glue, epoxy, tape or other suitable 
30 attachment means to ensure suitable contact between the sensor and the pipe 12. The 
sensors may alternatively be removable or permanently attached via known 
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mechanical techniques such as mechanical fastener, spring loaded, clamped, clam 
shell arrangement, strapping or other equivalents. Alternatively, the strain gages, 
including optical fibers and/or gratings, may be embedded in a composite pipe. If 
desired, for certain applications, the gratings may be detached from (or strain or 
5 acoustically isolated from) the pipe 12 if desired. 

Referring to Figs. 29,30, it is also within the scope of the present invention 
that any other strain sensing technique may be used to measure the variations in strain 
in the pipe, such as highly sensitive piezoelectric, electronic or electric, strain gages 
attached to or embedded in the pipe 12. Referring to Fig. 29, different known 

10 configurations of highly sensitive piezoelectric strain gages are shown and may 
comprise foil type gages. Referring to Fig. 30, an embodiment of the present 
invention is shown wherein pressure sensors 14-18 comprise strain gages 320. In this 
particular embodiment strain gages 320 are disposed about a predetermined portion of 
the circumference of pipe 12. The axial placement of and separation distance AXi, 

15 AX2 between the pressure sensors 14-18 are determined as described herein above. 

Referring to Figs. 31-33, instead of measuring the unsteady pressures P1-P3 on 
the exterior of the pipe 12, the invention will also work when the unsteady pressures 
are measured inside the pipe 12. In particular, the pressure sensors 14-18 that measure 
the pressures Pi,P2,P3 may be located anywhere within the pipe 12 and any technique 

20 may be used to measure the unsteady pressures inside the pipe 12. 

Referring to Figs. 34-36, the invention may also measure the speed of sound 
of a mixture flowing outside a pipe or tube 425. In that case, the tube 425 may be 
placed within the pipe 12 and the pressures P1-P3 measured at the outside of the tube 
425. Any technique may be used to measure the unsteady pressures P1-P3 outside the 

25 tube 425. Referring to Fig. 34, for example, the tube 425 may have the optical wraps 
302-306 wrapped around the tube 425 at each sensing location. Alternatively, any of 
the strain measurement or displacement, velocity or accelerometer sensors or 
techniques described herein may be used on the tube 425. Referring to Fig. 35, 
alternatively, the pressures P1-P3 may be measured using direct pressure 

30 measurement sensors or techniques described herein. Any other type of unsteady 
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pressure sensors 14-18 may be used to measure the unsteady pressures within the pipe 
12. 

Alternatively, referring to Fig. 36, hydrophones 430-434 may be used to sense 
the unsteady pressures within the pipe 12. In that case, the hydrophones 430-434 may 
5 be located in the tube 425 for ease of deployment or for other reasons. The 

hydrophones 430-434 may be fiber optic, electronic, piezoelectric or other types of 
hydrophones. If fiber optic hydrophones are used, the hydrophones 430-434 may be 
connected in series or parallel along the common optical fiber 300. 

The tube 425 may be made of any material that allows the unsteady pressure 
10 sensors to measure the pressures P,-P 3 and may be hollow, solid, or gas filled or fluid 
filled. One example of a dynamic pressure sensor is described in co-pending 
commonly-owned US Patent Application, Serial No. (Attorney Docket No. 712- 
2.40/CC-0067) entitled "Mandrel Wound Fiber Optic Pressure Sensor", filed June 4, 
1 999. Also, the end 422 of the tube 425 is closed and thus the flow path would be 
15 around the end 422 as indicated by lines 424. For oil and gas well applications, the 
tube 425 may be coiled tubing or equivalent deployment tool having the pressure 
sensors 14-18 for sensing P,-P 3 inside the tubing 425. 

Referring to Fig. 17, there is shown an embodiment of the present invention in 
an oil or gas well application, the sensing section 51 may be connected to or part of 
20 production tubing 502 (analogous to the pipe 12 in the test section 51) within a well 
500. The isolation sleeve 410 may be located over the sensors 14-18 as discussed 
hereinbefore and attached to the pipe 502 at the axial ends to protect the sensore 14-18 
(or fibers) from damage during deployment, use, or retrieval, and/or to help isolate the 
sensors from acoustic external pressure effects that may exist outside the pipe 502, 
and/or to help isolate ac pressures in the pipe 502 from ac pressures outside the pipe 
502. The sensors 14-18 are connected to a cable 506 which may comprise the optical 
fiber 300 (Fig. 22,23,27,28) and is connected to a transceiver/converter 510 located 
outside the well 500. 

When optical sensors are used, the transceiver/converter 510 may be used to 
receive and transmit optical signals 504 to the sensors 14-1 8 and provides output 
signals indicative of the pressure P,-P 3 at the sensors 14-18 on the lines 20-24, 
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respectively. Also, the transceiver/ converter 5 1 0 may be part of the Fluid Parameter 
Logic 60. The transceiver/converter 510 may be any device that performs the 
corresponding functions described herein. In particular, the transceiver/ converter 510 
together with the optical sensors described hereinbefore may use any type of optical 
5 grating-based measurement technique, e.g., scanning interferometric, scanning Fabry 
Perot, acousto-optic-tuned filter (AOTF), optical filter, time-of-flight, and may use 
WDM and/or TDM, etc., having sufficient sensitivity to measure the ac pressures 
within the pipe, such as that described in one or more of the following references: A. 
Kersey et aL, "Multiplexed fiber Bragg grating strain-sensor system with a Fabry- 

10 Perot wavelength filter", Opt. Letters, Vol 18, No. 16, Aug. 1993, US Patent No. 
5,493,390, issued Feb. 20, 1996 to Mauro Verasi, et al., US Patent No. 5,317,576, 
issued May 31, 1994, to Ball et al., US Patent No. 5,564,832, issued Oct. 15, 1996 to 
Ball et al., US Patent No. 5,513,913, issued May 7, 1996, to Ball et al, US Patent No. 
5,426,297, issued June 20, 1995, to Dunphy et al., US Patent No. 5,401,956, issued 

15 March 28, 1995 to Dunphy et al., US Patent No. 4,950,883, issued Aug. 21, 1990 to 
Glenn, US Patent No. 4,996,419, issued Feb. 26, 1991 to Morey all of which are 
incorporated by reference. Also, the pressure sensors described herein may operate 
using one or more of the techniques described in the aforementioned references. 

A plurality of the sensors 10 of the present invention may be connected to a 

20 common cable and multiplexed together using any known multiplexing technique. 

It should be understood that the present invention can be used to measure fluid 
volume fractions of a mixture of any number of fluids in which the speed of sound of 
the mixture a m j X is related to (or is substantially determined by), the volume fractions 
of two constituents of the mixture, e.g., oil/water, oil/gas, water/gas. The present 

25 invention can be used to measure the speed of sound of any mixture and can then be 
used in combination with other known quantities to derive phase content of mixtures 
with multiple (more than two) constituents. 

Further, the present invention can be used to measure any parameter (or 
characteristic) of any mixture of one or more fluids in which such parameter is related 

30 to the speed of sound of the mixture a m i X , e.g., fluid fraction, temperature, salinity, 
mineral content, sand particles, slugs, pipe properties, etc. or any other parameter of 



-31 - 



WO 00/00793 



PCT/US99/14589 



the mixture that is related to the speed of sound of the mixture. Accordingly, the 
logic 48 (Fig. 1) may convert a™ to such parameter(s). 

Further, the invention will work independent of the direction of the flow or the 
amount of flow of the fluid(s) in the pipe, and whether or not there is flow in the pipe. 
5 Also, independent of the location, characteristics and/or direction(s) of propagation of 
the source of the acoustic pressures. Also, instead of a pipe, any conduit or duct for 
carrying a fluid may be used if desired. 

Also, the signals on the lines 20,22,24 (Fig. 1) may be time signals 
Hi(t),H 2 (t)3 3 (t), where Hn(t) has the pressure signal Pn(t) as a component thereof, 
10 such that FFT[H,(t)] = G(fi>)P,(a>), FFT[H 2 (t)] = G(<o)P 2 (co), and the ratio 

H 2 (co)/H,(a)) - G(a>)P 2 (o)/G(cD)P,(©) = P 2 (a>)/P,(©), where G(a>) is a parameter 
which is inherent to each pressure signal and may vary with temperature, pressure, or 
time, such as calibration characteristics, e.g., drift, linearity, etc. 

Also, Instead of calculating the ratios P, 2 and P, 3 , equations similar to Eqs. 
15 9,10 may be derived by obtaining the ratios of any other two pairs of pressures, 

provided the system of equations Eq.5-7 are solved for B/A or A/B and the ratio of 
two pairs of pressures. Also, the equations shown herein may be manipulated 
differently to achieve the same result as that described herein. 

Still further, if, for a given application, the relationship between A and B (i.e., 
20 the relationship between the right and left travelling waves, or the reflection 

coefficient R) is known, or the value of A or B is known, or the value of A or B is 
zero, only two of the equations 5-7 are needed to determine the speed of sound. In that 
case, the speed of sound a mix can be measured using only two axially-spaced acoustic 
pressure sensors along the pipe. 
25 Further, while the invention has been described as using a frequency domain 

approach, a time domain approach may be used instead. In particular, the Eqs. 5,6,7 
may be written in the form of Eq. 1 in the time-domain giving time domain equations 
Pi(xi,t), P 2 (x 2 ,t),P 3 (x 3 ,t), and solved for the speed of sound a mix and eliminating the 
coefficients A,B using known time domain analytical and signal processing 
30 techniques (e.g., convolution). 
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Referring to Figs. 37-40, it should be understood that although the invention 
has been described hereinbefore as using the one dimensional acoustic wave equation 
evaluated at a series of different axial locations to determine the speed of sound, any 
known technique to determine the speed at which sound propagates along a spatial 
5 array of acoustic pressure measurements where the direction of the source(s) is (are) 
known may be used to determine the speed of sound in the mixture. The term acoustic 
signals as used herein, as is known, refers to substantially stochastic, time stationary 
signals, which have average (or RMS) statistical properties that do not significantly 
vary over a perdetermined period of time (i.e., non-transient ac signals). 

10 For example, the procedure for determining the one dimensional speed of 

sound a m ix within a fluid contained in a pipe using an array of unsteady pressure 
measurements is similar to a problem encountered in underwater acoustics (e.g., 
SONAR or Sound Navigation Ranging). In underwater acoustics, axial arrays of 
sensors are deployed to determine the bearing (or direction) of underwater noise 

15 sources. The process is referred to as "beam forming". In free space, i.e., in an 
unbounded media, such as the ocean, the speed at which a sound wave propagates 
along an axial array is dependent on both (1) the free-space speed of sound and (2) the 
incident angle of the sound wave on the axial array. 

Referring to Fig. 37, the apparent sound speed a x at which the wave 

20 propagates along the array is related to the angle or bearing (0 = 90 - y) of the source 
SI and the sound speed a in the media. For a SONAR application, as is known, the 
speed of sound is known and the apparent sound speed a x is measured, which allows 
the bearing to be determined by the relation: 0 = cos' 1 (a/a x ). 

Conversely, referring to Fig. 38, we have found that in a pipe 12 where the 

25 angle or bearing on the array of the incident sound is known, i.e., 0 = 0 deg, the speed 
of sound a of the fluid in the pipe 12 can be determined as follows. 

In particular, referring to Fig. 39, for a single distant source in two 
dimensional (2D) space, the pressure wave can be written as follows (such as is 
generally described in A. Dowling and J. Williams, "Sound and Sources of Sound", 

30 Ch4,pp79-81): 



-33- 



WO 00/00793 



PCT/US99/14589 



P(x,y 9 t)= Ae ia){t - xsinr * ,a - ycosr * ,a) Eq27 

Pressure as seen on the array at y=0 is: 

P(x,y = 0,t) = Ae ia{ '- xsinr > U) ^ 
P(x,t) = Ae- ik - x e i0 » Eq . 29 

5 where: k xr - (sin y x ) — . 

a 

A similar analysis may be done for a left travelling wave along the array from 
the source S2 as: 

P(x 9 t) = Be +ik « x e i ' 0 ' Eq.30 

where: k xl = (sin y 2 ) — . 

a 

10 For the situation where the sound is propagating along a pipe, then yi= y 2 = 90 

deg. and where a = amix which is the speed of sound of the fluid mixture in the pipe, 
then: 

7 7 ® 

k xr=Ki= Eq.31 

Thus, referring to Fig. 38, for a left and right travelling acoustic waves 
15 travelling in the pipe 12, the pressure equation becomes: 

P(x,t) = Ae- ik " x e imt + Be +ik * x e imt Eq .32 

which is the same as Eq. 1, and which may be used to determine the speed of sound 
by using the sensors described herein and solving the associated equations Eq. 5-7 
shown hereinbefore. The same result may also be shown from sources originating in 

20 three dimensional space using cylindrical or other coordinate systems. 

The data from the array of sensors may be processed in any domain, including 
the frequency/spatial domain (such as Eq. 4), the temporal/spatial domain (such as Eq. 
1), the temporal/wave-number domain or the wave-number/frequency (k-a>) domain. 
As such, any known array processing technique in any of these or other related 

25 domains may be used if desired. 
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For example, Eq. 5 can be represented in the k-co domain by taking the spatial 
Fourier transform of Eq. 5, resulting in the following k-co representation: 



P{k^) ^—Ipix^e^dx = A(co)S(k - — ) + B(a))s(k + — 
2x_i a \ a) 



Eq. 33 



where k is the wave number and 8 is the Dirac delta function, which shows a 
5 spatial/temporal mapping of the acoustic field in the k-co plane. 

Alternatively, instead of using the three equations Eq. 5-7, any technique 
known in the art for using a spatial (or phased) array of sensors to determine the 
direction of an acoustic source in three dimensional sound field with a known speed 
of sound (e.g., spatial array processing for SONAR arrays, RADAR (RAdio Detecting 
10 And Ranging) arrays or other arrays, beam forming, or other signal processing 

techniques), may be used to solve for the sound speed knowing the direction of travel 
of the acoustic waves, i.e., axially along the pipe. Some of such known techniques are 
described in the following references, which are incorporated herein by reference: H. 
Krim, M. Viberg, "Two Decades of Array Signal Processing Research - The 
15 Parametric Approach", IEEE Signal Processing Magazine, pp 67-94, R. Nielson, 
"Sonar Signal Processing", Ch. 2, pp51-59. 

Referring to Fig. 40, accordingly, the fluid parameter logic 60 may comprise 
spatial array processing logic 450 which receives the spatial array of acoustic pressure 
signals Pi(t), P2(t), P3(t) and performs the spatial array processing described herein to 
20 determine the speed of sound a m j X on the line 46. 

It should be understood that any of the features, characteristics, alternatives or 
modifications described regarding a particular embodiment herein may also be 
applied, used, or incorporated with any other embodiment described herein. 

Although the invention has been described and illustrated with respect to 
25 exemplary embodiments thereof, the foregoing and various other additions and 
omissions may be made therein and thereto without departing from the spirit and 
scope of the present invention. 



-35- 



WO 00/00793 



PCT/US99/14589 



Claims 

What is claimed is: 

1. An apparatus for measuring at least one parameter of a mixture of at least one 
5 fluid in a pipe, comprising: 

a spatial array of at least two pressure sensors, disposed at different axial 
locations along the pipe, and each measuring an acoustic pressure within the pipe at a 
corresponding axial location, each of said sensors providing an acoustic pressure 
signal indicative of the acoustic pressure within the pipe at said axial location of a 
10 corresponding one of said sensors; and 

a signal processor, responsive to said pressure signals, which provides a signal 
indicative of a speed of sound of the mixture in the pipe. 

2. The apparatus of claim 1 wherein said signal processor comprises logic which 
15 calculates a speed at which sound propagates along said spatial array. 

3. The apparatus of claim 1 wherein said signal processor comprises logic which 
calculates a frequency based signal for each of said acoustic pressure signals. 

20 4. The apparatus of claim 2 wherein said signal processor comprises logic which 
calculates a ratio of two of said frequency based signals. 

5. The apparatus of claim 1 comprising at least three of said sensors. 

25 6. The apparatus of claim 1 comprising three of said sensors and wherein said signal 
processor comprises logic which simultaneously solves the following equations for 
said speed of sound: 

P(x l9 t) = (Ae^ + Be + **)e im 
P(x 2y t) = (Ae + Be }?'*' 
30 P(x, 9 t) = (Ae^ +Be+ ik ' x >y° ,t 
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7. The apparatus of claim 1 wherein said signal processor calculates said speed of 
sound of said mixture using the following relation: 



o' ik r x \ 



+ Re 



+ Re 



3(g) 
p 3 (fi>) 



= o 



where i? = — = 



where £ r = 
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\a mix )\ + M x 



and &/ = 



CO 



10 where a m i X is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), 
and M x is the axial Mach number of the flow of the mixture within the pipe, where: 



15 



and where Vmix is the axial velocity of the mixture, and where Pi(co),P2(co),P 3 (co) are 
said frequency based signals for each of said acoustic pressure signals. 

8. The apparatus of claim 1 wherein said sensors are equally spaced, a Mach number 
of the mixture is small compared to one, and said signal processor calculates the 
speed of sound of the mixture using the following relation: 

0) 



a. = 



" 1 " 
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20 where P12 = Pi(co)/P2(co), P13 = Pi(co)/P3(co), 1 is the square root of-1, Ax is the axial 
spacing between sensors, where a mix is the speed of sound of the mixture in the pipe, 
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10 



15 



20 



a> is frequency (in rad/sec), and where Pi(co)J > 2 (o)) 3 P 3 ((o) are said frequency based 
signals for each of said acoustic pressure signals. 

9. The apparatus of claim 1 wherein said sensors are equally axially spaced, a Mach 
number of the mixture is small compared to one, and said signal processor calculates 
the speed of sound of the mixture using the following relation: 



Ax is the axial spacing between said sensors, and where Pi(co),P 2 (co),P3(q)) are said 
frequency based signals for each of said acoustic pressures signals. 

10. The apparatus of claim 1 wherein the signal processor comprises logic which 
calculates a fluid composition of the mixture in the pipe. 

11. The apparatus of claim 1 wherein said signal processor comprises logic which 
calculates the fluid composition of the mixture using the following relation: 



where ai,a 2 are known speeds of sound, p h p 2 are known densities, and hi,h 2 are 
volume fractions of the two respective fluids, and a mix is the speed of sound of the 
mixture. 

12. The apparatus of claim 1 wherein said speed of sound is substantially determined 
by two fluids within the mixture. 

13. The apparatus of claim 12 wherein said two fluids are: oil/water, oil/gas, or 
water/gas. 




where a^x is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), 
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14. The apparatus of claim 1 wherein said pressure sensors are fiber optic pressure 
sensors. 

5 15. The apparatus of claim 1 wherein at least one of said pressure sensors comprises 
a fiber optic Bragg grating-based pressure sensor. 

16. The apparatus of claim 1 wherein at least one of said pressure sensors measures a 
circumference-averaged pressure at said axial location of said sensor. 

10 

17. The apparatus of claim 1 wherein at least one of said pressure sensors measures 
pressure at more than one point around a circumference of the pipe said given axial 
location of said sensor. 

15 18. The apparatus of claim 1 wherein at least one of said pressure sensors measures 
strain on the pipe. 

19. A method for measuring at least one parameter of a mixture of at least one fluid 
in a pipe, said method comprising: 
20 measuring acoustic pressures within the pipe at at least two predetermined 

axial measurement locations along the pipe; and 

calculating a speed of sound of the mixture using said acoustic pressure 
measured at said axial measurement locations. 

25 20. The method of claim 19 wherein said calculating step comprises calculating a 
speed at which sound propagates along said axial measurement locations. 

21. The method of claim 19 wherein said calculating step comprises calculating a 
frequency based signals for said acoustic pressures. 

30 
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22. The method of claim 21 wherein said calculating step comprises calculating ; 
ratio of two of said frequency based signals. 



23. The method of claim 19 wherein said measuring step comprises measuring 
5 acoustic pressure at at least three axial measurement locations along the pipe. 

24. The method of claim 19 wherein said measuring step comprises measuring 
acoustic pressure at three axial measurement locations along the pipe and wherein 
said calculating step comprises simultaneously solving the following equations for the 

10 speed of sound: 



iat 



P(x x ,t) = (Ae~ ik ' x i + Be + *> X >}> 
P{x 2 ,t) = (Ae- ik ^ +Be + *' x >)e' 
P(x Jt t) = (Ae "*'*» + Be + ik,x > )?'' 



15 25. The method of claim 19 wherein said calculating step calculates said speed of 
sound of the mixture using the following relation: 



e -ik r x 3 +Re ik,x 3 



P.ico) 



= 0 



20 



where Rs — = 
A 



g-*r*, _ 


/}(*>)" 






P 2 (a>) 




P t (a>) 


e** - 




lP 2 (<o) 







where k r s 



<v 1 
N — 77~ md k s 



(O 



\ a mix)l-M x 



where a mix is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), 
and M x is the axial Mach number of the flow of the mixture within the pipe, where: 
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V . 



and where Vmix is the axial velocity of the mixture, and where Pi(cd),P 2 (q)),P3(o)) are 
said frequency based signals for each of said acoustic pressures. 

26. The method of claim 19 wherein said measurement locations are equally axially 
spaced, a Mach number of the mixture is small, and said calculating step calculates 
the speed of sound of the mixture using the following relation: 

0) 



' 1 " 




i log 





2P, 



13 



where P] 2 = Pi(co)/P2(co), P i3 = P](co)/P 3 (co), i is the square root of -1, Ax is the axial 
10 spacing between sensors, where a mix is the speed of sound of the mixture in the pipe, 
o is frequency (in rad/sec), and where Pi(co),P2(a>),P3(a)) are said frequency based 
signals for each of said acoustic pressures. 



15 



20 



27. The method of claim 19 wherein said measurement locations are equally axially 
spaced, a Mach number of the mixture is small compared to one, and said calculating 
step calculates the speed of sound of the mixture using the following relation: 



P 2 (co) 



a 



mix J 



where a m i X is the speed of sound of the mixture in the pipe, o is frequency (in rad/sec), 
Ax is the axial spacing between said measurement locations, and where 
Pi(o>),P2(co),P3(co) are said frequency based signals for each of said acoustic pressures. 



28. The method of claim 19 further comprising calculating a fluid composition of the 
mixture in the pipe. 
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29. The apparatus of claim 19 further comprising calculating a fluid composition of 
the mixture using the following relation: 



a mix = 



1 + 



\ °\ pi h \ 4 

where a,,a 2 are known speeds of sound, p„p 2 are known densities, and h,,h 2 « 
volume fractions of the two respective fluids, ^ is the speed of sound of the 



30. The method of claim 19 wherein the speed of sound is substantially determined 
by two fluids within the mixture. 

31. The method of claim 30 wherein said two fluids are: oil/water, oil/gas, or 
water/gas. 



32. The method of claim 19 wherein said measuring step is performed by fiber optic 
pressure sensors. 



33. The method of claim 19 wherein said measuring step is performed by fiber optic 
Bragg grating-based pressure sensors. 

34. The method of claim 19 wherein said measuring step measures a circumference- 
averaged pressure at said axial location of said sensor. 

35. The method of claim 19 wherein said measuring step measures pressure at more 
than one point around a circumference of the pipe at said axial location of said sensor. 
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